Condensate fraction of molecules for a spin mixture of ultracold fermionic atoms 
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The condensate fraction of molecules for ultracold Fermi gases is investigated for the magnetic 
field below the Feshbach resonant magnetic field. Assuming that there is no loss of particles and 
energy during the adiabatic magnetic-field sweep, a simple theory is used to interpret the measured 
condensate fraction in the experiments by JILA group (Phys. Rev. Lett. 92, 040403 (2004)) and 
MIT group (Phys. Rev. Lett. 92, 120403 (2004)). Our theory shows that the condensate fraction 
of molecules is dependent on the initial condition of the system and especially on the process of the 
magnetic-field sweep. 
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The evidence for Bose-Einstein condensates of di- 
atomic molecules has been finally given in several re- 
markable experiments 0, 0, which will obviously lead 
to intensive theoretical and experimental researches on 
the ultracold Fermi gases. The magnetic-field Feshbach 
resonance 0, Q plays an important role in most of the 
recent experiments on ultracold Fermi gases because it 
can change both the strength and sign of the scattering 
length a between fermionic atoms with different inter- 
nal freedom. On the side of strongly repulsive interac- 
tion (BEC side), there is molecule which is short-range 
fermionic pairs. On the side of strongly attractive inter- 
action (BCS side), one expects that there are fermionic 
pairs analogously to the electronic Cooper pairs in super- 
conductor. The Feshbach resonance has given us an im- 
portant opportunity to investigate the pairing phenom- 
ena in the ultracold Fermi gases especially the BCS-BEC 
crossover which has been discussed in a lot of theoreti- 
cal researches 0, Q S 0, EH E3, E3, E3, EM IS E3l 

Recently, the BCS-BEC crossover is inve stig ated by sev- 
eral experiments [3 E |24| such as the 
superfluid hydrodynamics p2l |23| and pairing gap 

The condensate fraction near a Feshbach resonance has 
been investigated experimentally in [LsL by adiabati- 
cally sweeping the magnetic field to different value, and 
the condensate fraction is then measured by a following 
fast decreasing of the magnetic-field to convert the dimers 
of fermionic atom pairs into bound molecules. Above the 
Feshbach resonant magnetic field, the condensate frac- 
tion observed by JILA group is investigated in ^3 
based on the model of a molecular BEC. On the Fesh- 
bach resonance, a simple theoretical model of a mixture 
of Fermi degenerate gas and dimeric gas p"?| is proposed 
to explain the high condensate fraction observed in [l9| . 
Below the Feshbach resonant magnetic field, the conden- 
sate fraction of molecules is investigated experimentally 
by JILA group m for 40 K and MIT group [HI for e Li 
with different magnetic-field sweep process. As far as 
we know, below the Feshbach resonant magnetic field, 
presently there is no theoretical interpretation for the 
condensate fraction of molecules measured in mm. in 
the experimental work of JILA group ^3 > the condensed 



atom pairs was found to be at most 15%, while the MIT 
group |l9j| observed a much higher fraction of condensed 
atom pairs below the Feshbach resonant magnetic field. 
The high condensate fraction observed by MIT group has 
been confirmed by other group such as the experiment by 
Bartenstein et al [2(j for 6 Li. Thus, there is an urgent 
need for a theoretical model to interpret this great differ- 
ence for different groups. In the present work, omitting 
the loss of particles and thus the loss of energy during 
the adiabatic magnetic- field sweep, a very simple theo- 
retical model is used to interpret the experimental results 

in both m for 4 ° K and - 19 ) for 6ll lt is found that the 
condensate fraction of molecules is dependent on the ini- 
tial condition of the system and especially on the process 
of the magnetic-field sweep. 

The magnetic-field Feshbach resonance 0, Q can 
change both the strength and sign of the atomic interac- 
tion. Near the Feshbach resonance, the scattering length 
is a (B) = dbg [1 — w/ (B — Bo)] with being the back- 
ground scattering length and B denoting the magnetic 
field with Feshbach resonant magnetic field Bo. Below 
the Feshbach resonance, the energy of the molecular state 
can be estimated well to be ~h 2 /ma 2 |2j| with m the 
mass of the fermionic atom. 

During an evaporative cooling process, it is well-known 
that there is a significant loss of energy and particles. For 
most of the recent experiments, however, the BCS-BEC 
crossover is investigated for the ultracold gases prepared 
at an initial magnetic field after different evaporative 
cooling process. After the evaporative cooling and dur- 
ing the following adiabatic magnetic-field sweep, the loss 
of particles can be omitted and thus we omit the loss of 
the energy of the system. In the present experiments on 
the BCS-BEC crossover, the temperature is much smaller 
than the Fermi temperature. Thus, in this work, we in- 
vestigate the ultracold gases at zero temperature. At 
zero temperature and during an adiabatic magnetic-field 
sweep process, omitting the loss of particles and thus 
the loss of energy, in the present work a simple theory 
is used to investigate the condensate fraction during the 
adiabatic magnetic-field sweep. 

Assume that N denotes the total number of fermionic 
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atoms in the absence of molecules and N m is the total 
number of molecules. Assume further that and 
are the number of fermionic atoms with different internal 
freedom. In the present experimental works on the BCS- 
BEC crossover, an equal mixture of fermionic atoms with 
different internal freedom is prepared in an optical trap. 
In this case, one has the following confinement condition: 



N = 27Vf + N, 



(1) 



To use the assumption that the loss of energy can be 
omitted, we calculate the overall energy of the system 
below and on the Feshbach resonant magnetic field. Be- 
low the Feshbach resonant magnetic field with \a\ /I < 1 
(J is the average distance between atoms), the chemical 
potential of the Fermi gas is given by 



gas in the unitarity limit can be given by the dimension- 
ality analysis: 



/4 es = (! + /?) 
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(6) 



where (3 is the correction due to the strong interaction 
between fermionic atoms with different internal freedom 
on resonance. The parameter has been inve stig ated 
by both experiments [20L I2H l30| and theories |26t l27l . 
On resonance, a quantum Monte Carlo calculation [27| 
gives (3 = —0.56. From Eq. JfjJ), the overall energy of the 
system on resonance is given by 



E res (N) = (1 + /3) 1/2 Wfj^-huho 



(7) 



(6^ 2 ) 2/3 ft 2 2/3 

nJ 

2m 



+ 7:9 + V F , (2) 



where the coupling constant gp — Ai:h 2 a/m, and n-f is 
the density distribution of the fermionic atoms in an in- 
ternal freedom. Vp = m(uj 2 x 2 + uj 2 y 2 4- lu 2 z 2 )/2 is the 
external potential of the fermionic atom. The overall 
energy Ef {N^N^ , B) of the interacting Fermi gases is 
given by 



E F (N h N u B) 



(6jV T ) 4/3 fuj h0 



+ g F n'dV, (3) 



1/3 

where uiho = (w z w s w z ) is the geometric average of the 
harmonic frequencies. 

At zero temperature, by using the Thomas-Fermi ap- 
proximation |3l|. the overall energy of N m molecules in 
the condensate is given by 



E m {N m ,B) 



Nr, 
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(4) 



where the energy of the molecular state is e 
—h 2 /ma 2 . In the above expression, 

\2/5 



length a m ss 0.6a 



/2 with the molecular scattering 
and the harmonic oscillator length 
d'£ = ^Jh/2muJho- The overall energy of the mixture 
gases of fermionic atoms and molecules is then 



Ef-t 



E F (N h N l ,B)+E m (N m ,B). 



(5) 



On resonance, the absolute value of the scattering 
length a is divergent and thus much larger than the 
average distance 7 between particles. In this situation, 
the ultracold gas can be described in the unitarity limit 
[26l I27I I2H I23 where there is a universal behavior for 
the system. By using the local density approximation, 
the general form of the chemical potential for the Fermi 



We turn to discuss the condensate fraction during the 
adiabatic magnetic-field sweep with different initial con- 
dition of the system. Firstly we consider the case that 
the initial system in thermal equilibrium is prepared on 
resonance, and then the magnetic field is lowered adi- 
abatically below the Feshbach resonant magnetic field 
Bo. In this situation, if the loss of the energy during 
the magnetic-field sweep is omitted, one has the follow- 
ing confinement condition on the condensate fraction 



E F (N h N l ,B)+E m (N m , B) = E res (N) . (8) 

Combining with Eq. one can get the condensate 

fraction of molecules x = 2N m /N at zero temperature 
for the magnetic field below Bq. 

In the experiment by JILA group 0] for 40 K, the scat- 
tering length is a = 174cto (1 — 7.8/ (B — Bo)) with ao be- 
ing the Bohr radius and Bo = 202. 1G being the Feshbach 
resonant magnetic field. In addition, according to the ex- 
periment in [l3|, N « 2 x 10 5 and uj x /2tt = ujy/2-K = 320 
Hz, ui z — uj x /79. In |l8j| . the system was initially pre- 
pared at a magnetic field 235.6 G far above the Feshbach 
resonance. The magnetic field was then swept across Bo 
adiabatically and especially there was a holding time of 
about 3 ms on resonance (see the inset of figure 2 in |18|). 
The adiabaticity of this sweep process can be verified by 
estimating the relaxation time 7 ~ \j\f2n\av with the 
total cross section a — 8iTa 2 and the mean value of the 
velocity of the thermal atoms v — ^SksT/ivm. A sim- 
ple calculation shows that for the regime of the magnetic 
field investigated in the relaxation time is smaller 
than 0.1 ms, while the sweep and holding time in [Tsj 
is larger than 1 ms. One should note that the conden- 
sate fraction was measured in |l8j | by lowering quickly 
(non-adiabatically) the magnetic field which results in a 
measured loss of 50% of the molecules . 

In 01, there was a holding time of about 3 ms on 
resonance. On the other hand, the divergent scattering 
length on resonance means that the relaxation time is 
extremely small and a small loss of particles on resonance 
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• experimental result for K by JILA group 
--0.56 
/3=-0.56 
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□ experimental result for Li by MIT group 
theoretical result of the c ondensate fraction 




FIG. 1: By using the experimental parameters of JILA group 
[Isj . shown is the condensate fraction x = 2N m /N for 40 K2 
molecules at zero temperature and below the Feshbach res- 
onant magnetic field Bo- The parameter j3 = —0.56 is the 
correction to the energy of the ultracold gases on resonance. 
The dotted line is the condensate fraction obtained from Eq. 
®, while the solid line is the condensate fraction obtained 
by considering the loss of molecules during the probing pro- 
cedure. The solid circles are the measured condensate fraction 
by JILA group. 

will change the overall energy of the system effectively. 
Thus the energy of the ultracold gases on resonance can 
be determined by E res (N). In this situation, Eqs. Q 
and JSJl can be used to calculate the condensate fraction 
x for different magnetic field B (< Bq). The dotted line 
in figure 1 shows the condensate fraction calculated from 
Eq. JHJ by using the above experimental data. We see 
that at least qualitatively the condensate fraction based 
on (3 = —0.56 agrees with the experimental result shown 
by the solid circle. The solid line shows the condensate 
fraction when 50% loss of the molecules is considered 
during the probing procedure (i.e. lowering the magnetic 
field quickly) pointed out in 18] . 

We now consider another case that the system is pre- 
pared initially below the Feshbach resonant magnetic 
field (i.e. on the side of molecular BEC). Assuming that 
Xini is the condensate fraction at an initial magnetic field 
Bi n i , the condensate fraction after lowering adiabatically 
the magnetic field to a value B is then determined by the 
following confinement condition when the loss of energy 
is omitted: 

E F (iV-j- = N i = N (1 - x ini ) A B ini ) + E m (x ini N/2, B mi ) 



= E F (iVf = N L = N (1 - x) /2, B) + E m (Nx/2, B). 

_ (9) 

Different from the experiment by JILA group [l|| , the 
experimental results illustrated by □ in figure 3 of the 



FIG. 2: The solid line shows the calculated condensate frac- 
tion of molecules at zero temperature for an experimental 
process of MIT group |Tg| where the ultracold gases with con- 
densate fraction 0.56 were initially prepared at the magnetic 
field 770 G below the Feshbach resonant magnetic field Bo- 
We see that the solid line based on the energy conservation 
given by Eq. agrees well with the experimental result 
shown by □ with a hold time of 100 ms. 



experimental paper by MIT group |li3| were observed 
for the system initially prepared at the magnetic field 
Bi n i = 770 G |l9j which is smaller than Bq. From the ex- 
perimental data, the condensate fraction at B ini = 770 G 
is Xi n i = 0.56. The condensate fraction of molecules be- 
low Bo was then investigated experimentally by lowering 
adiabatically the magnetic field. One can easily verify the 
adiabaticity of the magnetic-field sweep in |lj)] through 
the estimation of the relaxation time. The experimental 
parameters in [l]| are a = -1018a (1 + 325/ (B - B )) 
HH with B = 822 G and N = 2 x 10 6 , lu x /2tt = uj v /2it = 
115 x V25 Hz, u z /2ir = 1.1 x ^25 + 120 x B/1000 Hz. 

The solid line in figure 2 shows the theoretical result 
based on the above equation ©. We see that our the- 
ory agrees well with the experimental result of |19| illus- 
trated by □ in figure 2. Different from the analyses for 
the experiment by JILA group, we do not consider the 
loss of molecules during the probing procedure because 
in the experiment by MIT group for 6 Li, the molecules is 
much more stable comparing with the JILA group with 
40 K . We believe that both the experimental data of JILA 
group 0| and MIT group 0] are reliable. The large dif- 
ference of the condensate fraction in and is due 
to different experimental process. 

In the above theoretical investigation for the experi- 
ment by MIT group, the harmonic frequency lo z is de- 
pendent on the magnetic field. In fact, due to the in- 
dependent harmonic frequency u z , there is an effective 
work on the system by the magnetic field due to the con- 
tinuous changing of the ground state of the system. This 
effective work can be omitted because Aw z /w z < 4% 
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between the regime of magnetic field 770 G and 700 G 
investigated in the experiment. In figure 2, the agree- 
ment between theoretical result and experimental data 
shows clearly that the one can safely omit the effective 
work done by the magnetic field. 

For the condensate fraction discussed here, the term 
N m e m in the overall energy of the molecular gas 
E m {N m , B) plays a special role because it is negative 
and dependent strongly on the magnetic field. Below the 
Feshbach resonance and with the decreasing of the mag- 
netic field, the term e m decreases significantly. When 
the loss of energy is omitted, to maintain the energy con- 
servation, there should be a decreasing of the number of 
molecules when lowering adiabatically the magnetic field, 
which has been shown clearly in and ■ O ne should 
note that the overall energy of the ultracold gases is in- 
vestigated here for the case of weakly interacting gases 
with a/7 < 1 and the case of the divergent scattering 
length. Thus the theoretical result for the condensate 
fraction obtained here does not hold in the strongly in- 



teracting gases with a/7 > I, and this is the reason why 
we do not give the condensate fraction for a/7 > 1. 

In summary, the simple theory developed here is used 
to interpret different experiments by JILA group |18| and 
MIT group |l9j on the condensate fraction. Obviously, 
the condensate fraction discussed here is dependent on 
the magnetic-field sweep process and initial condition of 
an experiment. When the ultracold gases are prepared 
at the magnetic field below Bq, after different evapora- 
tive cooling one can get relatively arbitrary condensate 
fraction. One can check our theory further by investigat- 
ing the condensate fraction in experiments for different 
magnetic-field sweep process and different initial condi- 
tion of the system. It will be an interesting future work 
to investigate the condensate fraction of the fermionic 
atom pairs above the Feshbach resonant magnetic field 
based on the assumption that the loss of energy can be 
omitted. 

H. W. Xiong acknowledges very useful discussions with 
J. E. Thomas and C. Chin. 
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